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Abstract We propose a phantom crossing Dvali-Gabadadze-Porrati (DGP) model. 
In our model, the effective equation of state of the DGP gravity crosses the phantom 
divide line. We demonstrate crossing of the phantom divide does not occur within 
the framework of the original DGP model or the DGP model developed by Dvali and 
Turner. By extending their model, we construct a model that realizes crossing of the 
phantom divide. We find that the smaller the value of the new introduced parameter 
j3 is, the older epoch crossing of the phantom divide occurs in. Our model can account 
for late-time acceleration of the universe without dark energy. We investigate and show 
the property of Phantom crossing DGP model. 

Keywords modified gravity • Extra dimensions ■ equation of state • phantom 
crossing ■ cosmic acceleration • cosmological model 



1 Introduction 

Late time accelerated expansion of the universe was indicated by measurements of dis- 
tant Type la supernovae (SNe la) [TJ{9] . This was confirmed by observations of cosmic 
microwave background (CMB) anisotropies by the Wilkinson Microwave Anisotropy 
Probe (WMAP) 10, 11 , and the large-scale structure in the distribution of galaxies 
observed in the Sloan Digital Sky Survey (SDSS) [T5irT3] . 

It is not possible to account for this phenomenon within the framework of general 
relativity containing only matter. Therefore, a number of models containing "dark en- 
ergy" have been proposed as the mechanism for the acceleration. There are currently 
many dark energy models, including cosmological constant, scalar field, quintessence, 
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and phantom models |14H22) . However, dark energy, the nature of which remains un- 
known, has not been detected yet. The cosmological constant, which is the standard 
candidate for dark energy, cannot be explained by current particle physics due to its 
very small value, and it is plagued with fine-tuning problems and the coincidence prob- 
lem. 

An alternative method for explaining the current accelerated expansion of the uni- 
verse is to extend general relativity to more general theories on cosmological scales. 
Instead of adding an exotic component such as a cosmological constant to the right- 
hand side (i.e., the energy-momentum tensor) of Einstein's field equation, the left-hand 
side (i.e., the Einstein tensor, which is represented by pure geometry) can be modified. 
Typical models based on this modified gravity approach are f(R) models 23 - 25 and 
the Dvali-Gabadadze-Porrati (DGP) model [26H28] (for reviews, see [29]). 

In f(R) models, the scalar curvature R in the standard Einstein-Hilbert gravi- 
tational Lagrangian is replaced by a general function f(R). By adopting appropriate 
function phenomenologically, f(R) models can account for late-time acceleration with- 
out postulating dark energy. 

The DGP model is an extra dimension scenario. In this model, the universe is 
considered to be a brane; i.e., a four-dimensional (4D) hypersurface, embedded in a 
five-dimensional (5D) Minkowski bulk. On large scales, the late-time acceleration is 
driven by leakage of gravity from the 4D brane into 5D spacetime. Naturally, there is 
no need to introduce dark energy. On small scales, gravity is bound to the 4D brane 
and general relativity is recovered to a good approximation. 

According to various recent observational data including that of Type la supernovae 
30 33 , it is possible that the effective equation of state parameter w e ff, which is the 
ratio of the effective pressure p c g to the effective energy density p e g, evolves from being 
larger than —1 (non- phantom phase) to being less than —1 (phantom phase |16II34| ); 
namely, it has currently crossed —1 (the phantom divide). 

f(R) models that realize the crossing of the phantom divide have been studied 
35, 36 . On the other hand, in the original DGP model 26 - 28 and a phenomenological 
extension of the DGP model described by the modified Friedmann equation proposed 
by Dvali and Turner [37], the effective equation of state parameter never crosses the 
w g = — 1 lineQ 

In this paper, we develop the "Phantom Crossing DGP model" by further extend- 
ing the modified Friedmann equation by Dvali and Turner [37]. In our model, the 
effective equation of state parameter of DGP gravity crosses the phantom divide line, 
as indicated by recent observations. 

This paper is organized as follows. In the next section, we summarize the original 
DGP model, and check the behavior of the effective equation of state. In Section O 
we describe the modified Friedmann equation by Dvali and Turner [37], and we also 
demonstrate that the effective equation of state does not cross the iu e gp = — 1 line in 
this framework. In Section 2] we construct "the Phantom Crossing DGP model" by 
extending the modified Friedmann equation proposed by Dvali and Turner. We show 
that the effective equation of state parameter of our model crosses the phantom divide 
line, and investigate the properties of our model. Finally, a summary is given in Section 


1 Some models realize phantom crossing by adding an exotic component such as the scalar 
field to the DGP model | 38II39I . In these models, the equation of state of the additional 
component crosses the phantom divide, but the effective equation of state of DGP gravity 
itself does not. 
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2 The DGP model 

The DGP model [26] assumes that we live on a 4D brane embedded in a 5D Minkowski 
bulk. Matter is trapped on the 4D brane and only gravity experiences the 5D bulk. 
The action is 

S = lL M ^ / d5x V-9(5) R (5) 
1071 Jbulk 

+ T7T M U) / d ix J~9(i){R(i) + L m ), (1) 
1071 Jbrane 



where the subscripts (4) and (5) denote quantities on the brane and in the bulk, 
respectively. (M^) is the 5D (4D) Planck mass, and L m represents the matter 

Lagrangian confined on the brane. The transition from 4D gravity to 5D gravity is 
governed by a crossover scale r c . 

rc =2#- (2) 

On scales larger than r c , gravity appears 5D. On scales smaller than r c , gravity is 
effectively bound to the brane and 4D Newtonian dynamics is recovered to a good 
approximation. r c is the single parameter in this model. 

Assuming spatial homogeneity and isotropy, a Friedmann-like equation on the brane 
is obtained as [271128] 

2 87rG H 

= — p+t v c > (3) 

where p is the total cosmic fluid energy density on the brane. e — ±1 represents the two 
branches of the DGP model. The solution with e = +1 is known as the self-accelerating 
branch. In this branch, the expansion of the universe accelerates even without dark 
energy because the Hubble parameter approaches a constant, H — l/r c , at late times. 
On the other hand, e = — 1 corresponds to the normal branch. This branch cannot 
undergo acceleration without an additional dark energy component. Hence in what 
follows we consider the self-accelerating branch (e = +1) only. 

For the second term on the right-hand side of Eq. 0, which represents the effect 
of DGP gravity, the effective energy density is 

3 H (A\ 

Pr ^8^G- c ' (4) 



and the effective pressure is 



where H = dH/dt, the differential of the Hubble parameter with respect to the cos- 
mological time t. Using Eqs. and l[5). the effective equation of state parameter of 
DGP gravity is given by 

w rc = ^. (6) 

pr c 



Fig. \T\ shows the behavior of the effective equation of state of DGP gravity «v c 
versus the redshift z for r c Ho = 1.04, 1.43, and 2.00. Assuming that the total cosmic 
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fluid energy density p of Eq. @ contains matter and radiation, from Eq. (J7J, these 
values of r c H$ correspond to J? TO .o = 0.04, 0.30, and 0.50, respectively. 

Q mfi = 1 - C2 rfi - (rcHoy 1 , (7) 

where f} m is the normalized energy density of matter and Q r is the radiation on the 
brane; i.e., Q m = (8nG/3H 2 )p m and Q r = (8nG/3H 2 )p r . (p m«a p r oc a ^). The 
subscripts designate the present value. 

The effective equation of state of DGP w rc can also be exactly expressed in terms 
of the energy densities of matter and radiation, fl m and Q r [401141] . 

W Tr = . (8) 

In realistic ranges of the energy density, (2 m > and O r > 0, the value of the effective 
equation of state cannot be less than or equal to —1. That is, the effective equation of 
state never crosses the phantom divide line in the original DGP model. 



3 DGP model extended by Dvali and Turner 



Dvali and Turner [37] phenomenologically extended the Friedmann-like equation (Eq. 
((3}) of the DGP model. This model interpolates between the original DGP model and 
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the pure ylCDM model with an additional parameter a. The modified Friedmann-like 
equation is |37] 

2 87rG H a , . 

For a = 1, this agrees with the original DGP Friedmann-like equation, while a — 
leads to an expansion history identical to that of ylCDM cosmology. 

Differentiating both sides of Eq. ([9| with respect to the cosmological time t, we 
obtain the following differential equation. 

2H = -8*G(p+P)+ ^"y^ ( io) 

where a dot indicates the derivative respect to the cosmological time. The quantity P 
is the total cosmic fluid pressure on the brane. 

For the second term on the right-hand side of Eq. @, which represents the effect 
of DGP gravity, the effective energy density is 

Pa = ^-pr 2 _ a , 11) 



and from Eq. (|10|l the effective pressure is 

(12) 



(rcH) 



From Eqs. (|11[) and (|12|) . the effective equation of state parameter of the DGP model 
extended by Dvali and Turner is given by 

W a = — . (13) 
pa 

Fig. [2] shows a plot of the behavior of the effective equation of state of the DGP 
model by Dvali and Turner w a versus the redshift z for a — 1.50, 1.00, 0.50, 0.00, —0.50, 
and —1.00 (assuming i? m ,o = 0.30). 

In general, for equation of state w, the energy density p varies as a~ 3< -- 1+w \ This 
leads to the following proportional relation. 

p a oc a . (14) 

At the same time, From Eq. (JTTJ), we find p a oc H a . In the radiation-dominated epoch, 
from the 
relation. 



from the proportional relation on Hubble parameter H oc a 2 , we obtain the following 



pa oc a 2q . (15) 

As compared the right-hand side of Eq. (|14|l to that of Eq. I|15|l . during the earlier 
radiation-dominated epoch [z ^> 10 ), the effective equation of state can also be rep- 
resented with a |37| . 

2a ,,„, 
w a = -1 + —■ (16) 

In the same way, during the matter-dominated epoch (10 2 > z3>l), from the propor- 

3 

tional relation on Hubble parameter H oc a 2 , 

Wa = -l+2. (17) 
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Fig. 2 Effective equation of state w a , vs. redshift z of the DGP model extended by Dvali and 
Turner for a = 1.50, 1.00, 0.50, 0.00, -0.50, and -1.00 (top to bottom) assuming Q mfi = 0.30. 



At the present time, w a is a stationary value close to —1. 

From these results, in the case of a < 1.0, the effective equation of state w a < — 3 
in all era, even in the radiation-dominated epoch. That is to say, the component of the 
DGP gravity works as the driving force of the accelerated expansion of the universe 
in all epochs. On the other hand, for a ^. 1.0, there is era when the effective equation 
of state becomes w a ^ — 3. That is, the DGP gravity does not drive the accelerated 
expansion in all epochs. 

The case of a — 1.0 corresponds to the original DGP model described in the 
previous section. Thus, in the original DGP model, the effective equation of state 
Wr c — — j in the radiation-dominated epoch. And after the radiation-dominated epoch, 
becomes w rc < —3- In other words, the DGP gravity acts as the driving force of the 
accelerated expansion just after the radiation-dominated epoch. 

However, when a is positive, the effective equation of state w a will exceed —1 at 
all times. For negative a, w a is always less than —1. In the case of a = 0, w a is — 1 
constantly. Based on this analysis, crossing of the phantom divide does not occur in 
the DGP model extended by Dvali and Turner. 
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4 Phantom crossing DGP model 



We propose the "Phantom Crossing DGP model" that extends the modified Friedmann 
equation (Eq. (J2|) proposed by Dvali and Turner. Our model can realize crossing of 
the phantom divide line for the effective equation of state of the DGP gravity. 

As mentioned in the previous section, the effective equation of state parameter of 
the DGP model by Dvali and Turner w a , takes the value of over — 1 for positive a, 
and it is less than —1 for negative a. When a — 0, w a becomes — 1. On the basis of 
these results, we consider a model in which a varies being positive to being negative. 
To keep the model as simple as possible, we make the following assumption, 

a = P-a, (18) 

where a is the scale factor (normalized such that the present day value is unity). The 
quantity f3 is a constant parameter. In the period when the scale factor a is less than 
the parameter f3 (a > 0), the effective equation of state exceeds —1. At the point when 
the scale factor a equals /?, (a = 0), the equation of state's value will be —1. In the 
period when the scale factor a exceeds the parameter f3 (a < 0), the equation of state 
will be less than — 1. In this way, crossing of the phantom divide is realized in our 
model. 

Replacing a by /3 — a in Eq. @, the Friedmann-like equation in our model is given 

by 

H =—P+ rc 2-(e-a) - ( 19 ) 

Differentiating both sides of Eq. (|19p with respect to the cosmological time t, the 
following differential equation is obtained. 

a r(, i P\ i ii - a )H~aH In (r c H) 



For the second term on the right-hand side of Eq. (|19|l representing the effect of DGP 
gravity, the effective energy density is 



H 



P-a 



and from Eq. (|20[) . the effective pressure is 



1 

"8^G 



(/3 - a)H - aH\n(r c H) #< 3 ~ a 



(21) 



(22) 



(rctf) 2 -^-") rc 2-(/3-a) 
Using Eqs. (|21|) and Q22p . the effective equation of state of our model is given by 

yp = ^r- (23) 



Fig. [3]shows a plot of the effective equation of state of our model wp versus the red- 
shift z (see also Fig. [4] which shows an enlarged view of this diagram). Our model is an 
extension of the DGP model and realizes crossing of the phantom divide. The effective 
equation of state wp of models for /3 = 0.50,0.25, and 0.10 (assuming ft m fi = 0.30) 
crosses the phantom divide line when the redshift z ~ 0.2,0.8, and 1.6, respectively. 
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We find that the smaller the parameter B is, the older epoch crossing of the phantom 
divide occurs in. 8 is not necessarily equal to the scale factor at the time of crossing 
the phantom divide, even though Eq. ()18|) is assumed. In the Eq. (|19p . the value of 
B — a that is the power index of H varies with respect to time. As the power index 
of the differential equation changes with time, furthermore, in parallel, the differential 
equation is solved with respect to time. Hence, the time lag occurs, the scale factor at 
the time of crossing the phantom divide is more than the value of 8. 

In a way similar to the derivation of Eq. (|16p . we represent the effective equation of 
state of Phantom Crossing DGP model with B. In the radiation-dominated epoch, the 
scale factor a is taken to be in comparison with the value of B. Therefore, as a ~ B 
in Eq. (JTSJ), during the radiation-dominated epoch (z S> 10 4 ), the effective equation of 
state is approximately 

28 

ifl/jsa-l+y. (24) 

That is, in the case of B < 1.0, the effective equation of state wp < — jjj in all era, 
including the radiation-dominated epoch. On the other hand, for B 2; 1.0, there is era 
when the effective equation of state becomes wp ^ — jj. 

Fig. [S] shows the effective energy density pa and absolute value of the effective 
pressure IP^I (note that Pp < 0) of our model for (B, f2 m ,o) = (0.50,0.30) versus 
the redshift z, normalized such that the effective energy density is unity at the time 
of phantom crossing. It shows that the absolute value of the effective pressure \Pp\ 
exceeds the effective energy density pp at the time of crossing of the phantom divide. 
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Fig. 5 The red (solid) and pink (dotted) lines respectively represent the effective energy 
density pa and absolute value of the effective pressure | Pa \ (note that Pa < 0) of our model 
vs. redshift z, for {fi, Q m>0 ) = (0.50,0.30). 
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The recent observational data for Type la supernovae 5 show that crossing of the 
phantom divide line occurs at a redshift z ~ 0.2 30 32 01n our model, for f3 = 0.50 
(when f2 m fi = 0.30), crossing of the phantom divide occurs at z ~ 0.2. 

In a proposed model in which the phantom divide is crossed at z ~ 0.2, (/?, f2 m ,o) 
— (0.50, 0.30), we investigate and show the property of Phantom crossing DGP model. 




z 

Fig. 6 The distance modulus fj, relative to that of a constant expansion cosmology fi c , vs. 
the redshift z. Models and parameters are (from top to bottom): (1) ylCDM model, O m = 
0.30; (2) Phantom Crossing DGP model, = 0.50, Q m .o = 0.30; (3) DGP model by Dvali'and 
Turner, a = 0.50, J? m , = 0.30; (4) Original DGP model, i? mj0 = 0.30. 



Fig. [6] shows the distance modulus fj, relative to that of a constant expansion cos- 
mology fi c , versus the redshift z. That is, when fi/fic is positive, cosmic expansion is 
accelerating. The distance modulus is defined by 

p{z) = 5 log 10 D L + 42.38 - 5 log 10 h, (25) 

where Dl is the Hubble free luminosity distance given by 

D L = (l + z)£^^, (26) 

2 This is a model-dependent value. We will investigate in detail the allowed parameter region 
based on recent observational data in future work. 
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z 

Fig. 7 The normalized energy density of radiation i? r , matter Q m , and DGP gravity Qdgp> 
vs. the redshift z in the Phantom Crossing DGP model with the proposed parameter (/3, f2 m ,o) 
= (0.50,0.30). 



h being the Hubble constant Hq in units of 100 km s^ 1 Mpc -1 . We adopt h = 0.72 [42]. 
In Fig. [5] Models and parameters are (from top to bottom): (1) ylCDM model, J? m .o 
= 0.30; (2) Phantom Crossing DGP model, f) = 0.50, Q mja = 0.30; (3) DGP model 
by Dvali and Turner, a = 0.50, J? m , = 0.30; (4) Original DGP model, J? m , = 0.30. 
Phantom Crossing DGP model can realize late-time acceleration of the universe very 
similar to that for ylCDM model, without dark energy. 

Fig. [7] shows the normalized energy density of radiation i? r , matter Q m , and DGP 
gravity f2oGP versus the redshift z in the Phantom Crossing DGP model with the 
proposed parameter (j3,Q m fi) = (0.50,0.30). Where, Odgp = {8ttG/3H 2 )p^. pp is 
the effective energy density of DGP gravity defined by Eq. (|2ip . We find that the 
universe is DGP gravity-dominated near recent epochs. Therefore, In the Phantom 
Crossing DGP model, the late-time acceleration is driven by the effect of DGP gravity. 



Fig. |8]shows the effective equation of state ui e // versus the redshift z. Models and 
parameters are same as Fig. [|J] Only our Phantom Crossing DGP model can realize 
crossing of the phantom divide line at z ~ 0.2 as indicated by recent observations. 
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Fig. 8 The effective equation of state w e ff, vs. the redshift z. Models and parameters are 
same as Fig. [6] 



5 Summary 

— We confirmed that the effective equation of state does not cross the phantom divide 
line in the original DGP model. 

— We also demonstrated that crossing of the phantom divide does not occur in the 
DGP model by Dvali and Turner. 

— We constructed the Phantom Crossing DGP model. This model realizes crossing 
of the phantom divide. We found that the smaller the value of the new introduced 
parameter f3 is, the older epoch crossing of the phantom divide occurs in. Our model 
can realize late-time acceleration of the universe very similar to that of /1CDM 
model, without dark energy, due to the effect of DGP gravity. In the proposed 
model ((/3, fi m ,o) = (0.50, 0.30)), crossing of the phantom divide occurs at z ~ 0.2 
as indicated by recent observations. 
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